NASA Contractor Report 181928 
ICASE Report No. 89-47 


ICASE 


NEAR-PLANAR TS WAVES AND LONGITUDINAL 
VORTICES IN CHANNEL FLOW: NONLINEAR 
INTERACTION AND FOCUSING 


P. Hall 
F. T. Smith 

•, 1 1 

u ncl is 
\j? 3 7 u S j 

Contract No. NAS 1-1 8605 
September 1989 

Institute for Computer Applications in Science and Engineering 
NASA Langley Research Center 
Hampton, Virginia 23665-5225 

Operated by the Universities Space Research Association 


t, jH I \! AL V j -1 l T l -> 


L' ITU 
NO ML I N cAH 


AN 


i-ie po r * 


( I C A j r ) 


A,. T r WA V ' 

C H'N ' t ‘ 1 w * 

; r u'l. f N r rin-1 

C<-(A 01 A 


0J/0? 


NASA 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton, Virginia 23665-5225 




NEAR-PLANAR TS WAVES AND LONGITUDINAL 
VORTICES IN CHANNEL FLOW: NONLINEAR 
INTERACTION AND FOCUSING 

P. Hall 1 


Department of Mathematics 
University Exeter 
EX4 4QE, U.K. 


F. T. Smith 

Department of Mathematics 
University College London 
Gower Street, London WC1E 6BT, U.K. 


ABSTRACT 

The nonlinear interaction between planar or near-planar Tollmien-Schlichting waves and 
longitudinal vortices, induced or input, is considered theoretically for channel flows at high 
Reynolds numbers. Several kinds of nonlinear interaction, dependent on the input amplitudes 
and wavenumbers or on previously occurring interactions, are found and are inter-related. 
The first, Type a, is studied the most here and it usually produces spanwise focusing of both 
the wave find the vortex motion, within a finite scaled time, along with enhancement of both 
their amplitudes. This then points to the nonlinear interaction Type b where new interactive 
effects come into force to drive the wave and the vortex nonlinearly. Types c,d correspond 
to still higher amplitudes, with c being related to b, while d is connected with a larger-scale 
interaction e studied in an allied paper. Both c,d are subsets of the full three-dimensional 
triple-deck-lie interaction, f. The strongest nonlinear interactions are those of d,e,f since they 
alter the mean-flow profile substantially, i.e., by an 0(1) relative amount. All the types of 
nonlinear interaction however can result in the formation of focussed responses in the sense 
of spanwise concentrations and/or amplifications of vorticity and wave amplitude. 

1 This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NAS1-18605 while the first and second author was in residence at the Institute for Computer 
Applications in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 


1 




1. INTRODUCTION 


A fairly common feature in experiments and large-scale computations on channel-flow 
transition, as well as boundary-layer transition, with or without wall curvature, is the ap- 
pearance quite early on of significant streamwise streaks following the span wise concentration 
of streamwise vorticity. These streaks are possibly associated with the creation of Lambda 
vortices. Channel-flow experiments are described by Nishioka et al. (1975, 1980, 1985), 
Kozlov and Ramazanov (1982, 1983, 1984), Ramazanov (1985), among others, and related 
computations are by Kleiser et al. (1982, 1985), Orszag et al. (1980, 1983), Biringen et 
al. (1984), Zhou and Wang (1984), Zang and Hussaini (1985), Gilbert and Kleiser (1986), 
Singer et al. (1986), among others, under various disturbance conditions. It appears how- 
ever that no rational or other theory has been put forward yet to account for the focusing 
of streamwise vorticity and the eventual formation of streaks. The present theoretical work 
is aimed at providing such an account, as well as enlarging the understanding of Tollmien- 
Schlichting/longitudinal-vortex interaction. 

The theory is for an incompressible fluid in unsteady three-dimensional (3D) motion at 
high Reynolds numbers, starting with an initial input which is anear-2D Tollmien-Schlichting 
(TS) wave. A slight 3D “warping” present in the TS wave provokes a slight vortex motion as 
part of the mean-flow correction. This then poses a number of physical and theoretical ques- 
tions. The first question to be addressed is how nonlinear interaction can take place between 
the evolution of the TS wave and that of the induced vortex flow, i.e., what structures are 
involved? The second question concerns whether the interaction is self-sustaining or instead 
peters out, and, third, how dependent is it on the initial input? Fourth, can theory predict 
the change of scale and spanwise focusing that seem central for the eventual description of 
streak formation? 

These questions are considered in the following sections. In sections 2,3 it is found that 
the first occurrence (labelled Type a), or one of the first kinds, of such nonlinear interaction 
can arise for rather tiny amplitudes of the input TS wave 3 This tinyness is due essentially 
to the quite wide difference in the typical time scales for the viscous-inviscid TS wave and 
the induced vortex, which is initially viscous-controlled. The scales are derived from an 
order-of-magnitude argument in section 2, although they could be deduced instead from 
adaptation of the reasoning in Hall and Smith (1984, 1987). Next, the interaction turns out 
to be self-sustaining (section 4) in the sense that the slight warping of the 2DTS wave forces 
a very fast growth of the vortex motion which in turn reinforces the warping, and so on. At 
first this occurs in a form of linear secondary 3D instability which is different as a rule from 

Specifically, the TS streamwise velocity perturbation involved is of order Re -4 / 7 , suggesting a threshold 
of about 1% at Reynolds numbers Re in the range 1000-5000. 
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the suggested destabilization of Squire modes in that the secondary instability is found to be 
much more violent than previously suggested (although the subsequent interactions below do 
re-establish the importance of Squire- mode destabilization). After that, nonlinearity soon 
takes control and the indications, again from section 4, are that the nonlinear TS-vortex 
interaction terminates within a finite scaled time with the formation of concentrated zones 
or tongues of enhanced streamwise vorticity and TS amplitude (as well as phase variation), 
at particular spanwise locations. These vorticity tongues are felt to presage the initial de- 
velopment of a typical streak. A change of scale therefore takes place then, bringing in new 
physical properties associated with higher amplitudes inside the tongue, and this defines a 
second Type b of nonlinear interaction, as discussed in section 5. Further Types c,d,e are 
identified in section 6, along with the connections between the Types a-e. Throughout, the 
nonlinear behavior exhibits a sensitive dependence on the initial conditions; moreover, less 
civilized transitions than the experimental, computational and theoretical ones concerned 
above can by-pass any of the present regimes and enter later, higher-amplitude, stages di- 
rectly. Further comments on these and other aspects are presented in section 7. 

The velocities u,u, w and the corresponding cartesian coordinates x,y,z are based on 
nondimensionalization with respect to w#, one-quarter of the maximum speed of the undis- 
turbed plane Poiseuille flow, and to the channel width respectively. The undisturbed 
flow has u = u(y) = y — y 2 (with v, w zero) between the walls y = 0,1, properties at which 
are denoted by superscripts =p in turn. The pressure is written Pd^^dP and the Reynolds 
number is Re = upaojuD^ where po^D are the density and kinematic viscosity of the fluid. 
Time is written ant/un, and the subscripts r, i stand respectively for the real and imaginary 
parts of a quantity. The theory is concentrated nominally near the lower branch of the neu- 
tral curve for large Re: it happens that the Reynolds numbers of concern in the transition 
experiments are also large, and previous theoretical studies indicate that the present kind of 
approach can have application at relatively low, even subcritical, Reynolds numbers. 

2. THE SCALES, AND THE CORE FLOW, FOR TYPE A INTERACTIONS 

Only slight warping of a 2D TS wave in the incoming 2D plane Poiseuille flow is necessary 
to provoke a nonlinear interaction between the wave and the longitudinal vortices created 
by the warping. The scales involved arrange themselves as follows for the first, Type a, 
interaction to be considered here: see also Figure 1. First of all, the typical time or length 
scales of the primary TS wave and a longitudinal vortex are very different, the time scales 
being 0(e“ 3 ) and 0(e -7 ) in turn, where e = Re -1 / 7 is small. The former size comes straight 
from properties of the Orr-Sommerfeld equation at large Re and helps fix the thin viscous 
critical-layer/ wall-layer width, while the induced longitudinal vortex in the core responds 
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initially at least with an unsteady-viscous force balance, so that d t ~ Re -1 . This relative 
discrepancy of order £ 4 in the time scales or £ 6 in the lengths implies that a crucial amplitude 
of the 2D TS wave is of order £ 3 , in relative terms, since the critical nonlinear interplay 
is expected to be due to amplitude-squared inertial effects acting on the x-scale 0(£ -1 ) 
characteristic of a 2D TS wave, and that determines the TS disturbance amplitude in the 
core (0 < y < 1) and the wall layers (y ~ 0, 1) examined below. A similar reasoning based on 
physical balances also points to the slow z-scale of interest, namely z ~ e -3 , for the spanwise 
warping in the TS wave to be significant, as opposed to the main x-scale of 0(e _1 ). Likewise, 
the induced vortex amplitude follows as 0(£ 4 , in the velocity it, because of its production 
by means of amplitude-squared effects from the TS disturbance, in the core of the flow. So 
multiple scaling is called for, in the form 

d x — *■ edx, d z — y e 3 dzi d t — ► e 3 dr + e 7 df (2.1) 

where x = e 1 X, z = e~ 3 Z, t = e~ 3 T, t = e~ 7 T, and so on, along with the core- and 
wall-layer expansions set out just below and in the next section. A noteworthy feature in 
both cases is the fairly tiny size, in both the TS and the vortex disturbances, at which 
nonlinear interaction first enters; this tinyness seems to stem originally from their disparity 
in characteristic time scales and it seems to tie in with some experimental observations. 

In the core of the undisturbed motion, the flow solutions for the velocity and pressure 
take the form 

U = u(y) + £ 4 U (3) + £ V 1 ) + £ 7 li (e) + . . . + £ 9 U (/) + . . . 

V = e G V^ -f £ 8 lA e) + . . . + £ W V^ + . . . + £ 4 V 3 ) + . . . 

w = e s w ( - 3 ' 1 + £ 9 u> ( 1 ) + £ n u;( e ) + . . . + £ 13 iw^) -f . . . 

P = -Ke 6 X + e 7 pW + e 9 p( e ) + . . . + e n p^ + £ 12 p (3) + . . . 

where the superscript (3) refers to the induced vortex flow which is dependent on T, Z , y but 
independent of the fast scales X,T, while (1) refers to the nonlinear TS wave (dependent 
on X,T,T,Z,y) and (e), (/) denote higher-order TS effects. Also, the constant k 0 = 2 from 
the basic flow. Substitution of (2.1), (2.2) into the continuity and Navier-Stokes equations 
then leads to the following successive balances controlling the TS wave and the longitudinal 
vortex in the core. 

The TS wave at leading order satisfies, for continuity and the x — , y — , z — momentum 
balances in turn, 

Ujf(i) + tyi) = 0, uu X (i) + i/V = 0, 

«%(*) = -Pj/i), uw x( ,) = -p z( ,), 


(2.2a) 
(2.2 b) 

(2.2 c) 

(2.2 d) 
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(2.3a, b) 
(2.3 c, d) 



where the prime denotes d/dy. The solution here is neutral at this level of working, so that 

{u ( V (I) ,«> (1) ,P (1) } = {Av!,-iaAu, p (1 >, u>W}E + c.c. (2.4) 

where p (1) , w (1) are independent of X, T, the function A is independent of X , y, T and repre- 
sents a negative displacement effect, 


E = exp^aX - ITT)] (2.5) 

with the wavenumber a and the frequency 9 both being real, and (2.3a-c) are satisfied 
provided that 

P (1) = ~ <* A f u 2 dy, u? (1) = -p z(i )/iau. (2.6a, b) 

J o 

Here P {1) ~(X ,T,T, Z , . . .) represents the lower-wall pressure, and the scaled pressure differ- 
ence between the walls, 

pd)+ _ p(D- = _ a 2 AI [ wher e /= 1/30], (2.7) 

is due to the curvature in the disturbed core flow and is a vital part of the nearly 2D TS 
wave. Also, c.c. denotes the complex conjugate of the preceding expression. 

Next, the TS corrections . . . are governed by 


u x( «) T Vy( C ) — 0, u T ( i) + uu X ( C ) -f v^u' = — Px(i)-, (2.8a, b) 

u TO) + uv X ( t ) = -p y (e), wji i) + uw X ( C ) = -p Z ('), (2.8c, d) 

thus bringing in the unsteadiness and streamwise pressure gradient of the TS wave. A 
decomposition as in (2.4), 


{u 


(e) „(.) U ;W,pW} = (u(e) )i )(e) > (e)^(e)] £; + cc 


therefore gives the results 


(2.9) 


« (e) = ^ - u Jp {l) u~ 2 dy - p^/u, (2.10a) 

= iau J p^u~ 2 dy -f i9A — iaA^u, (2.106) 

P y («) = 2a9.Au + a 2 u 2 {-A^ + J p^u~ 2 dy} (2.10c) 

where — A^ is an extra displacement effect. 

Then the extra TS contributions u^\ ... are the first ones to feel the influence of the 
unknown induced vortex flow [u^, . . .] directly, through the governing equations 


u x(f ) + ty/) = 0, u T(c) + uu x( s) + u (3) u X (,) + v U) u > + u (1) tx y (3) = -p x( e) (2.11a, b ) 
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*>r(«) + UV X U) + W ( 3 > ^( 1 ) = -p y u),WTX*) + uw xU ) + u {3) w X ( q + v ( 1 ) wf J = -Pzu) (2.11c, d) 

and the nonlinear inertial interaction between (1) and (3) terms. We note that any terms 
between in ( 2 . 2 a) and so on do not affect the present vortex/TS interaction. With 

the decomposition = u^’E + c.c., etc., as in (2.4), (2.9), the solution here gives 

= —iat{A^u + Au^ + Ar^(i/)}, (2.12a) 

pW = pG)~ — a 2 {A ^ I u 2 dy + 2A f uu 3 dy} + A*y^\y) (2.126) 

Jo Jo 

where the functions are standard functions in the linear-stability expansion for a 

purely 2 D TS wave and need not concern us here, — is the displacement correction at 
this level and gives the corresponding lower-wall pressure contribution. Hence the 

small pressure difference induced across the channel at this order is given by 

/>(/)+ _ pi!)- = - a 2 {i (/) 7 + 2 AJ) + AK (2.13) 

were J = fo uu^dy and K = /J 7 ^dy. The result (2.13), as expected, is a slight distortion 
of the earlier pressure jump (2.7), due most significantly to the distortion from the vortex 
motion. 

Lastly here, the longitudinal-vortex flow etc., is controlled on the slower time scale 
by the equations 

ttys) + w Z ( 3 ) = 0, (2.14a) 

Uj*3)+ < + t! (1) *Vi) > + u(3) w' = « w (3), (2.146) 

< u (1) a x( i) -f > +G 0 «« (3) = ~p y ( 3 ), (2.14c) 

ra^T < u {l) w xW + t> ( 1 ) uyi) >= -p z w + w y ^)- (2.14c?) 

Here <> refers to the E° component of the enclosed nonlinear terms. Other components 
oc E , E 2 (T c.c.) are also produced but these do not influence the dynamics substantially. 
Again, Go is a scaled Gortler number whose effects could be included in the subsequent 
analysis although henceforth we take the case of zero Go to focus on flat- surface transition. 
So elimination of p^ and manipulation of (2.14) lead to the governing equations 

(d 2 y - df)dlv (3) = 4 a 2 u u'd z {\A\d z \A\} (2.15a) 

(d 2 y - df)u {3) = a'u (3) (2.156) 

of the longitudinal vortex motion in the core, subject to the boundary conditions 

u ( 3 ) _ y(3) _ Q y v^ = 0 at y = 0, 1 (2.15c) 
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for no slip at the channel walls, in view of (2.14a). There is an implicit assumption here that 
the induced vortices have little wall-layer structure/modification, and this is borne out by 
the balances acting in the wall layer as described next. 


3. THE WALL LAYERS IN TYPE A INTERACTIONS 


Viscous wall layers which are in effect critical layers are required near y — 0, 1 to reduce 
the TS flow, e.g., in (2.4), to rest at the walls. The previous scalings in the core suggest the 


expansions 

u = e 2 \ y + eW + e 5 C/ (1) + e 6 U {3) + e 7 t/ (e) + . . . + e 9 U {}) + ..., (3.1a) 

v = e 8 V (1) + £ 10 V (e) + • • • + £ 12 V (/) + • • • + e 15 V< 3 > + . . . (3.16) 

w = e 7 W {1) + e 9 W {e) + e 10 W (3) + . . . + £ n W U) + . . . (3.1c) 

p = -k 0 e e X + £ 7 P (1) + e 9 P {e) + ...+ £ n P (/) + £ 12 P (3) + . . . (3. Id) 


in the lower wall layer near y — 0, where y = e 2 y with y of 0(1), and similarly in the viscous 
upper wall layer of width 0(£ 2 ) adjoining y — 1. Here A = 1, p = — 1 from the basic flow, 
while U^ x \ U^ e \ are the induced TS contributions and = A3 _y is the induced-vortex 
streamwise shear flow from the core motion, so that A 3 - = A 3 -(T, Z) = d y u^ at y = 0-f. 
The vortex velocity components V^ 3 \ are less simple, being forced by TS-amplitude- 
squared effects in the wall layer, as in Hall and Smith (1987), although they match with the 
behavior oc y 2 ,y (plus logarithmic behavior) respectively at large y required by the core-flow 
results. Further, the y-variation in all the pressure terms shown is found to be negligible 
here, leaving i^ 1 ) — p( 1 )~E+ c.c., etc., independent of y. As before, the terms between 
the (e),(/) contributions in (3.1a) and so on, including second harmonics ex P 2 , have no 
significant impact on the vortex/TS interaction. 

The governing equations for U^ l \ for and for are basically the only ones required 
now. The first are unsteady, viscous, quasi-2D ones of a TS wave, 


V x(') + fy( j ) — 0? Ujii) + AyC/^(i) -f V^A — — P^o) + £^$0) (3.2a, b) 


subject to 

U {1) = V (1) = 0 at y = 0, U {1) -> XA as y -> oo, (3.2c) 

for no slip and for matching with the core solution. Hence U ^ = U^E+ c.c., etc., where 
from (3.2a-c) we have, with Ai denoting the Airy function, 


U$d = BAi{§, BAi ' 0 = iaPW A“§, 
A A = iap( l ) A~* k / Ai' 0 , 


(3.3a) 
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where £ = A *(y - 0 /A a), A = Am, the subscript zero denotes evaluation at £ = £ 0 = 
— Asfi/Ao and k = Ai(q)dq. The corresponding upper wall yields the relation 

-XA = iaP {1)+ A->K/Ai ' 0 , (3.36) 

on the other hand, due to the opposing displacement effect there, and so coupling (3.3a,b) 
with the core property (2.7) we obtain the dominant eigenrelation 

As K a 2 I = 2X 2 Ai' 0 (3.3c) 

fixing a, fi as both real. Computation gives the values a ~ ( 1 .00)(2/ /) 7 , Q ~ (2.298)as. 
Next, the relatively small spanwise velocity induced, satisfies 

Wf > A A yW ( ^ = -P ( z 1] + Wg (3.4) 

with VK (1) = 0 at y = 0 and Vb (1) -» 0 as y — ► oo. So kK (1) = W^E+ c.c. where 1T (1) is 
given by 

VT (1) = A"tp^ 1)_ L(0 (3.5) 

and the function L(£) is defined by L" -£L = 1, L 0 = L(oo) = 0, and is expressible in terms 
of integrals involving the Airy function. 

The extra TS contribution is then forced mainly by the slow variation of the primary 
TS wave, both spanwise in Z and temporally in T, and by the additional shear from the 
induced vortex motion. Other contributions are present from the correction U^ e \ V ^ in the 
TS wave but they are less significant. Thus the controlling equations are 

U xU ) A + W z(1) = 0, (3.6a) 

U TU) A A XyUxU) + yfU x(c) A t/< 3) t/ x(1) + V {1) U x a) A V (1) t/f 1 ' 

+2 V^ny A V (S) X = -P { X J) A (3.66) 

subject to the constraints 

ljU) — yU) — 0 at y = 0, — > A + A 3 - A + T^(0)A as y — * oo. (3.6c) 

The influences of the slow variation and the vortex shear, here and in the upper wall layer, 
must therefore satisfy a compatibility relation. This can be worked out by solving the system 
(3.6a-c) for its ^-components, then conducting a similar analysis on the forced contribution 
in the upper viscous layer, and relating the two via the cross-channel jump condition (2.13) 
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on the induced wall pressures. After some manipulation, and using the earlier solutions, we 
obtain the amplitude equation 


2( Aiq/ k)(£o T Aiq! k) As A t — — 2ot 4 A 3 J A A 2\ 2 (Ai 0 / A*c) A^z T ttj A 
— ('4*o/ K )(A 3 + + A 3 -)*a[l — |{^o 2 (Aio/i4? 0 ) — 1 + ioAiol «}]A 


(3.7) 


governing the effect of the induced vortex motion (through the terms J,A 3 ±) on the TS 
amplitude A(T,Z). In (3.7) the constant ttj is the standard linear-stability contribution for 
a 2D TS wave, determining whether the flow is sub- or super-critical, while the unknown 
integral property J(T, Z) and induced shears A 3 ±(7 1 , Z) are given by 

J = t uu^dy, A 3 - = d yU W{T, Z, 0), A 3+ = -d y u^\T, Z, 1). (3.8) 

Jo 

The equation (3.7) can also be checked against an appropriate expansion of the eigenrelation 
(3.3c), in fact. 

In summary, then, the nonlinear interaction between the warped 2D TS disturbance and 
the longitudinal vortex motion is controlled by the partial/ordinary differential-cum-integral 
system consisting of (2.15a-c), (3.7) along with the definitions of J, A 3 ± in (3.8). 


4. NONLINEAR TYPE A SOLUTION PROPERTIES; SPECIAL CASES; SEC- 
ONDARY 3D INSTABILITY; FORMATION OF FOCUSSING 


The system governing the nonlinear interaction may be written in the form: 



{d 2 y - df)d 2 y v (z) = 2 u v!d\p. 

(4.1a) 


(d 2 — dj)u^ — uv^\ 

(4.16) 

subject to 

u (3) _ ^(3) _ fl yV ( 3 ) = 0 at y = 0, 1 

(4.1c) 

and 

ckfA = [c 4 + c 2 J + c 3 (A^ T A 3 )]A T c^d^A, 

(4.14) 


J(T, Z) = f'uu^dy; A±(T,Z) = =p^« (3) at y = 1,0; p=\A\ 2 . 
Jo 

(4.1c) 


Here A(T,Z) = a A, while the constants c 2 ,c 3 ,c 4 follow from a computation for the coeffi- 
cients in (3.7) as 


c 2 = 131.5 - 175.2*', c 3 = -1.826 - 0.0348*', c 4 = -0.0655 + 0.0873*' (4.2) 
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and the constant Cx(<x Xx) has positive or negative real part depending on whether the basic 
flow u = y{ 1 - y) is super- or sub-critical. As might be expected, there are many features 
in common between the nonlinear system (4.1) and Hall and Smith’s (1987), we note, an 
exception being the term oc c 4 in (4. Id) which turns out to have a substantial impact. Again, 
there is also a connection with Bennett, et al.’s (1988), Hall and Smith’s (1989), and Smith’s 
(1988b) work on high-amplitude vortex effects when extreme Z-variations are considered. 
The main solution properties of interest are discussed in (a) - (g) below. 

(a) First, for an incident TS wave that is only relatively slightly warped we have \d z \ 
being small, say d z — with a small and Z typically 0(1). Initially then, for example, 
the vortex flow induced is small, 0(cr 2 ), and so A grows as an undistorted TS wave like 
exp(cxT) times an arbitrary function of Z, if Cj > 0, and the interaction is weakened. The 
0(<r 2 ) vortex motion then amplifies like exp(2cxT), however, being forced by the p zz term 
in (4.1a), and this drives a corresponding exp(3ciT) growth in the 0(cr 2 ) correction of the 
TS wave, via the interaction through J, Xf in (4. Id). Hence a fuller interaction comes back 

A - 

into play at increased times T of order —Ina, at which stage |A| is larger, of order a , but 
the vortex amplitude is also increased, to 0(1). Since the time variation remains typically 
0(1) the full interactive system (4.1a-e) is retrieved at that stage except that the c a A zz 
term in (4. Id) is absent due to the slower Z-variation here. This less-warped case therefore 
forms a connection with Hall and Smith’s (1987), and indeed it is found (see also the later 
computations and Figure 6) that the interactive system with C 4 in (4.2) replaced by zero 
yields nonlinear exponential growth of the TS and vortex amplitudes as T increases, in line 


with the findings of the last-named paper. Specifically, we have here 

[u (3) ,u (3) ,p,arg A] ~ e sf [v,u,^, <£] + ..., as T -> 00 , (4.3a) 

and the nonlinear growth factor s is determined by the reduced problem obtained from 
substitution into (4.1) with c 4 zero, namely 

(d 2 - s)dlv = 2 u vfdlp, (4.36) 

(dy — s)u = uv, (4.3c) 

u' = v = d y v — 0 at y = 0, 1, (4.3d) 

c 2r f uudy + c 3r {d y h(0, Z) - d y u( 1, Z)} = 0. (4.3c) 

Jo 


Here (4.3e) follows from the dominant terms of order exp(3sT/2) that appear in (4. Id) 
when (4.3a) holds. The profile p(Z) remains arbitrary, although we observe that (4.3a-e) 
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represents an exact solution of (4.1) when C\ = C\ = 0, in which event p(Z) is fixed by the 
initial distribution. The eigenvalue s is one of those determined by Hall and Smith (1987). 

(b) Next, suppose that C\ is small, corresponding to a state which is closer to neutral 
than the above, e.g., if c 4 is positive this occurs at a lower Reynolds number, for a prescribed 
wavenumber. Then the less-warped case initially has 

A = a 0 (f, Z) + a 2 Ai{T,T, Z) + ... (4.4a) 

where df — > df + a 2 df and 

Ayf + a 0 f = [ci + c 2 J + c 3 (A 3 + A 3 )]a 0 + QOqzz (4.46) 

from (4. Id), whereas (4.1a,b) lead to 

(d 2 y - df)d 2 y v {3) = 2 u up Q zz > ~ df)u {3) = uv {3) . (4.4c, d) 

Here 

= (4.4e) 

d = cr 2 Ci, and p 0 = |a 0 | 2 . Apart from transients, the correction term A\ grows exponentially 
with time [or a 0 does] and so at a later stage we have effectively the ordering 

A = 0(1), dz = &dz, df = a 2 df , (4.5) 


along with (4.4e). Thus the governing equations now reduce to 

d 2 y v {3) = 2 u up??, d 2 y u (3) = uv {3 \ (4.6a, 6) 

Af = [£i + c,J + c 3 (A 3 + A 3 )]A + c 4 A^zi (4.6c) 

since the time scale is relatively slow here. The boundary conditions on (4.6a,b) are that 
{|( 3 ), t)( 3 ), 5yi)(3) vanish at y = 0, 1 and the solution is straightforward, giving 

6 (3 > = {2 y 7 - 7 y 6 + 7r/ 5 - 3y 3 + y 2 }/420, (4.7a) 


u( 3 ) = { — y 10 /45 + y 9 /9 — 3y 8 /16 + y 7 / 1 2 + y 6 /10 
-y 5 / 8 + y 4 / 24 - l/720}/210 


(4.76) 


explicitly; hence the vortex properties J, A 3 can be evaluated and substituted into (4.6c). 
The resulting complex amplitude equation for A(T, Z) is then 



[ci + 61 pzz\A + c 4 Azzi 


(4.8a) 
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where 


(4.86) 


P=\M* 

61 = (-0.486 + 3.905i)10~ 5 . (4.9) 

We address this simplified nonlinear version, (4.8a,b), now because to some extent it mirrors 
the main properties of the full system (4.1). It is vaguely like certain well-studied equations 
such as the cubic Schrodinger equation, incidentally, although in general no firmer connection 
seems available. 

One particular solution is simply that for 2D flow of course where A = A 0 (f) is inde- 
pendent of Z, and there (4.8a) yields the exponential form A 0 = exp(c 1 T) in line with the 
definition of c x - To examine the stability of this 2D solution to small 3D distortions we 
use the polar form A = rexp(iO) with ( r,6 ) ( T,Z ) real, convert (4.8) to two real coupled 
nonlinear equations for r,0, 

r f ~ (^ir T 6i T pzz) r + c 4r( r zz ~ r ^z) — c 4i( r @z)z/ r \P = r2 ]i (4.10a) 

rOf = {hi + hiPzzY + c 4.( r zz - re \) + c 4r {r 2 0z)z/r (4.106) 

and then expand, with 8 small, 

M) = { 7 oe^ f A + hit] + 6{n,<h) + 0(6 2 ). (4.11a) 

The expression in curly brackets is the 2D solution, with 7o (> O),0 o prescribed constants, 
and of course p = p 0 + 8p x - f 0(6 2 ) where p 0 = 7 ^exp(2c lr T’),/> 1 = 2 7o rj exp (c lr f)- The 
linearized system resulting from substitution of (4.11a) into (4.10) yields solutions 

(ri,</q) oc exp(i(3Z)exp(Qf) (4.116) 

for small c Jr or, at large times, for c lr negative, and a similar result holds for c Xr positive. Here 
/? is the particular spanwise wavenumber of the 3D distortion and the effective growth-rate 
Q is found to be given by 

Q = hr ~ /? 2 (6 lr po + C 4r ) ± /5{6 IrPlfi 2 '- Abupohr ~ 2^ 2 b u p 0 c Ai - /3 2 c 2 4l }i. (4.11c) 

The dependence of Q r on (3 and the 2D wave amplitude />o is presented in Figure 2. Several 
points emerge here but the main ones are the following. First, the secondary 3D instability 
is very pronounced, especially for increased fi and/or p 0 where maxQr increases like /3 2 p 0 
approximately. Thus both the effects of three-dimensionality and of 2D amplitude growth 
act in concert to destabilize, accentuating the secondary distortion or warping. The next 
point is that the secondary 3D instability here is quite different as a rule from Herbert’s 
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(1984) suggested form and is associated more with the 3D growth-rate behavior at large 
Reynolds numbers, through the c 4r term above, than with the induced-vortex/Squire-mode 
response through bip 0 , c.f. Section 6 below. The latter effect dominates only if c 4 is negligible 
as in item (a) earlier corresponding to reduced three-dimensionality in the initial input, this 
giving a loose connection with Squire-mode destabilization. Although there is therefore some 
dependence on the initial conditions of course, in general the enhancement of any warping is 
due mainly to the c 4r contribution. Third, this can occur even while the flow is subcritical. 
Fourth, the increase in temporal growth rate with increasing is so fast that for an initial- 
value problem in general, where (4.11b) is replaced by a spectral decomposition in Z for 
instance and in effect integration with respect to /? becomes appropriate, the 3D distortion 
becomes singular within a finite time. This is similar to the break-up phenomenon associated 
with the backward heat-conduction equation. In the present case the break-up takes place 
with a time scale oc Q -1 or /3 -2 , typically, i.e., focusing occurs with \Z — Zs\ oc (7s — 7)2 at 
a finite time 7s and location Z$- The local behavior there is given by a similarity form which 
is analogous to that in (4.13) below and is somewhat dependent on the initial conditions. If 
these or the boundary conditions are unrestricted then the linearized initial-value problem 
can become ill-posed in the sense that effectively 7s tends to zero. This cumulative instability 
for large wavenumbers $ and/or large amplitudes p 0 is due in the present regime to & lr ,c 4r 
being negative but it permeates the analysis and computations in the other regimes studied 
as well. 

Another simple particular solution should also be mentioned, namely the single mode 
A = A 0 (7) exp(i/3Z) with 0 fixed, giving an extension of the previous one and yielding to 
similar analysis. Here in addition p = |A 0 | 2 = p(T) and so A 0 (T) satisfies the linear equation 

A' 0 = (C! - c 4 /3 2 )A 0 . (4.12) 

This further illustrates the strong instability possible since, as c 4r is negative, the solution 
oc exp((ci — c 4 /3 2 )7) grows exponentially for (— c 4r )/? 2 + Ci r positive and the growth rate 
increases indefinitely with increasing spanwise wavenumber /?. Limit-cycle solutions may 
also be sought but they are probably, although not necessarily, irrelevant in view of the more 
general behavior proposed next. 

The aspects above suggest that the most common response of the system (4.8a,b) is likely 
to be a nonlinear break-up within a finite time, T\ say. The orders of magnitude involved 
indicate that as 7 — ► 7\— the singular behavior has 

A ~ y- 2 a(r?), p ~ ^ p{rj ) and i] = |c 4 |2(Z - Z x )j(J x - 7)2 (4.13a - c) 

Oi Oj 
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focussed around the position Z = Z\. There, from (4.8a,b), the functions a(r]),p(ri) satisfy 


\r)a = p"aexp(i<j>i) + a"exp(?> 2 ), P - |«| 2 , 


(4.14a, 6) 


subject to a,/3 being bounded at large |r/|, to match with the flow away from Z « Zj, and 
with <^i = arg(6i), <^ 2 = arg(c 4 ) so that 7r/2 < </>i j2 < tt. The localized behavior here produces 
a shock-like jump from a(— oo) = ao to a(oo) = 5i, i.e., at Z = Zi±. For example, if the 
jump is small then (4.14a,b) yield the result 


/ rj 

ex p(-Nij 2 )drj 

-OO 


(4.15a) 


where 


N = [~{poC\ + c 2 )ci ± {(/3 0 C? - 5iS 2 ) 2 - S 2 }*] 


(4.156) 


/[4cj { 1 + 2po cos(</>i — ^2)}] 

and c n = cos <j) n <0, s n = sin <f> n > 0, n = 1,2. As each root N has positive real part the 
solution (4.15a) is acceptable and the proposed jump behavior is obtained from 7/ = —00 to 
77 = -|-oo. As regards other acceptable forms, numerical solutions of (4.14a,b) for a range 
of jump conditions were derived by use of a finite-difference procedure and these are shown 
in Figure 3. They agree well with (4.15) for small jumps. The description (4.13) - (4.14) 
of a local nonlinear break-up also gains weight from the fact that, at least for zero cj, it 
provides an exact similarity solution of the system (4.8a,b). The occurrence of the break-up 
here provides the means for the less-warped interaction to focus and acquire the full form 
(4.1a-e), leading on to the subsequent response studied in (f), (g) below. 


(c) The single-mode solution mentioned in (b) above may be extended to the full nonlinear 
system (4.1a-e), since if A — A 0 (T) exp(ifiZ) then A 0 is governed by (4.12) again provided 
Ci is replaced by ci. So the comments immediately following (4.12) apply also to the general 
case. 


(d) Likewise, secondary 3D instability can be established for the full system by extension 
of the approach in (4.10) - (4.11). For the full system (4.1) we use polars in the form 
A = Rexp(i<f>) to replace (4. Id) by the two real equations 


Rf — \c\ T + c 2r J + C3 r (A3 -f A 3 )]/? + C^ r [RzZ ~ R^z] 
—c 4i (R 2 <f)z)z/ R, 

R<pf = [ci, + c 2 , J -f c3i(A 3 + A 3 )]/? + c 4 i[Rzz — R<t> 2 z ] 
+C4 r{R 2 4>z)zl R 


(4.16a) 


(4.166) 
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with p = R? now, and then the secondary 3D instability of the 2D TS solution R = 
r 0 exp(c lr r), <j) = <j) Q -\- c u T [where To,(Ao are constants, To > 0] is examined by perturbing 
it, for small 6, 


(R, <f>) = {r 0 exp(ci r r), 4> 0 + c u t } + s(R l , (f> x ) + o(<5 2 ), 

= 6(uj 3 *,t4 3 ') +0(^ 2 ). 


(4.17a) 


The resulting linearized system for R x , (f>\ is found to yield both fast growth in the form 


R\ fa exp[7exp(2ci r T/3)] (7 = constant) 


(4.176) 


at large times, as might be expected, and fast growth for increasing spanwise wavenumbers 
I 3. The latter growth is analogous to that in item (b) above, causing a finite-time breakdown 
of the small 3D distortion in (4.17a) for an initial-value problem and making the subsequent 
computational work very sensitive. 


(e) An extreme case of interest opposite to those in (a), (b) arises for relatively large 
amplitudes and short spanwise lengths \Z\, such that Z is 0(T) say with T << 1. Then if 
|A| is 0(r _2 ),/9 is 0(r _1 ) and the balances in (4.1a-e) suggest that typically the time scale T 
is fast, 0(r 2 ), with v ^ ~ T~ 4 and id 3 ) ~ T -2 ~ A 3 ~ J all being large in magnitude. This 
has the effect of rendering the core flow of (4.1a,b) inviscid, apart from two wall layers of 
width 0(T) in y, which influence the values of A3 , and of replacing c x in (4. Id) by zero. The 
resulting simplified system then ties in closely with the analysis in (g) below. 

(f) Computational studies of the full nonlinear system (4.1a-e) (with (4.2)) have been 
made, applying both spectral and finite-difference approaches, for various starting conditions 
imposed at time T = 0 and for various values of the constant c x . IN the spectral treatment 
Fourier series in Z of period L were substituted into (4.1a-e) to yield a set of coupled 
nonlinear equations in y,T for the Fourier components. The set was truncated at a large 
number (A/) of terms and marched forward in time steps AT with a predictor-corrector 
method of second-order accuracy in y,T, for grid sizes Ay in y. A range of values for the 
parameters L, M , AT, Ay was also examined, and Figure 4 shows a typical result. Similarly, 
Figure 5 presents a typical computation from the finite-difference procedure, where time- 
stepping was applied as before but central differences were taken in both y and Z (step 
size A Z). The results from both methods tend to confirm the analytical features that we 
have described previously, with the properties of shortest spanwise scale growing the fastest 
temporally. This leads on to the account in (g) below. 

The reduced system implied in (a) above, i.e., (4.1a-e), (4.2) but with c 4 replaced by 
zero, has also been studied computationally and a representative result is given in Figure 6, 
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obtained from the spectral approach. The secondary 3D instability analogous to (4.11a-c) is 
much reduced now and the cascade of energy into higher spanwise components mentioned 
in the last paragraph is not observed; instead the nonlinear amplitude growth computed 
appears to be consistent with the description in (a) above, that is, the TS wave and the 
vortex flow amplify together exponentially. 

(g) The terminal behavior of the full system (4.1a-e), (4.2) is considered now. Many types 
of algebraic or exponential nonlinear breakdown suggest themselves and have been tried by 
us, 4 and the same goes for additional special cases, but the favorite description, guided by 
the features (a) -(f) above, is as follows. 

We propose that the nonlinear interaction breaks up by producing a concentrated “tongue” 
formation within a finite time T = T s : see Figure 7. The orders of magnitude suggest then 
the singular behavior 

i?=|i|~(T s -T)- 1 fi(5), p ~ (T. - T)~ 2 p(z), (4.18a, b) 

u< 3) ~ (T s -T)~ l u(y,z), ~ (T, -T)~ 2 v(y,z), (4.18c, d) 

<j> ~ Gln(T s -T) + (f>(z) (4. 18c) 

to leading order as T — > T s ~, with the dependence on y, z as shown, and with the tongue 
being focussed within a distance 0(T S — T) 2 of the spanwise location Z = Z 3 say, so that 

Z-Z S = (T,-T)*z (4.18/) 

where z is 0(1). The algebraic break-up (4.18a-f) can be generalized readily, but the par- 
ticular scales shown above appear to be the favored ones: c.f. (4.13). Also here, the vortex 
response (4.18c,d) applies in the bulk (i) (0 < y < 1) of the flow, where the vortex becomes 
mainly inviscid, whereas new viscous wall layers (ii) of width oc ( T a — T) 2 come into play 
near y = 0, 1 in which u^ 3 \v^ reduce in magnitude to the levels ( T s — T ) 2 , (T 5 — T) 2 
in turn. The vortex skin-friction contributions A3 are determined by the properties in those 
wall layers (ii) but the vortex’s coanda effect J is dominated by the contribution from the 
bulk (i), and so viscous-inviscid interaction maintains itself during the tongue development. 
The governing equations (4.1a,b) for the vortex motion become 

(—2 — zdz/2)dlv ” 2 u ufizz , (4.19a) 

(-1 -zd- z /2)u = u'v, (4.196) 

4 For example, the nonlinear exponential growth in (4.3) can be seen to apply in principle to the full 
system also, this growth then providing a match with our related work, Bennett, et al. (1988), Hall and 
Smith (1989), Smith (1988b) on the long-scale vortex/wave interactions. 
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(4.19c) 

(4.19c?) 


subject to u = v 


= 0 at y = 0, 1, while the TS part satisfies [with p = R 2 ] 


R + zR- z j 2 — [c 2 T J + C3 r ( + Ag )]R + Cir{Ra — R<j)\) ~ C4i(R 2 fa) | / R, 

R{z(t>i/2 - G) = [c2, J + c3,(Ag + A 3 )}R + C4 T (R 2 (j)z)z/ R + c^Ru — R(j>\) 
and, again to leading order, J ~ (T, — T) -1 J(z), A 3 ~ ( T s — T) _1 A 3 (z) with 


J = 



(4.19e) 


The skin friction terms A 3 , however, are affected by the wall layers (ii) which are required 
to reduce d y v to zero at t/ = 0,1. A solution of the terminal problem (4.19) may then be 
sought in the separated-variables manner: 


M] = [V a (y)V b (z), U a (y)U b (z)}. (4.20a) 

Thus (4.19a,b) split into the ordinary differential equations 

V” = 2 ku<n', (~2V b - zVj /2) = k~ l p" (4.206) 

U a = £xfv a , (~U b - zUi/2) = rV fc (4.20c) 

where U a , V a are to vanish at y = 0, 1, the prime denote differentiation with respect to y or z 
as appropriate and k,l are constants. Equations (4.19c,d) remain intact here provided that 
(J,A 3 ) is replaced by 

(J a , A^ )U b (z) with J a = f uU a dy. (4.20d) 

Jo 

The solution for V a is 

V a = k{[ u 2 dy ~ V (4.21a) 

from which U a follows in (4.20c), giving J a in (4.20d) equal to — H/4200. In contrast, V b , U b 
are given imphcitly by 

V b = —2k~ l z~ 4 [ Z z 3 p"{z)dz, U b = -2 r l z~ A [ Z z 3 V b {z)dz, (4.216, c) 

^0 Jo 

the origin of integration being zero to avoid an irregularity in general. Hence the terminal 
response may be reduced in principle to the solution of a pair of integro- differential equations 
for 72, 4> (using p — R 2 ) as functions of z in the range — oo < z < oo, by substitution of 
(4.20d) - (4.21c) into (4.19c,d). Local expressions valid for small \z\ in the center of the 
tongue can be written down and allow the solution to be symmetric there. In the other 
extreme, at the edges of the tongue as z — ► ±oo, the vortex and TS amplitudes tend to zero, 
specifically in the form 

R ~ RqZ 2 R\z poz 4 + p\z 6 + . . . , (4.22a, b ) 
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(4.22c) 
(4.22 d, e ) 


<j> ~ 2G£n\z\ -f- 2 + . . . , 

Vb ~ QqZ ~ 4 + a\Z ~ 6 + . . . , Ub ~ 2 + 4 + . . • , 

consistent with the controlling equations above. Here Ro,R\, Po, P\i$u a o, a ub 0 , are con- 
stants. The complete behavior in (4.22) then matches with the less pronounced flow outside 
the tongue where the amplitudes |A|, |u^|, |u^| and the phase <j) are all 0(1) at 0(1) distances 
from Z = Z,. 

Thus the streak pattern of (4.18a-f) appears to provide a self-consistent terminal response 
to the nonlinear interaction. The pattern is in fact another nominally exact solution of the 
full system (4.1a-e), (4.2) if the TS linear-growth term c 4 is neglected. A particular type of 
this terminal form we note has the TS amplitude dominating, equivalent to |«/|, |Af | being 
relatively small in (4.19c,d), in which case we find the solution 

Rexp(i<j)) = d x [2c A + zexp(z 2 /4c 4 ) f exp(-£ 2 /4c 4 )dz] (4.22) 

J d>2 

where d x ,d 2 are constants. This yields an acceptable tongue response, in agreement with 
(4.21), since c 4r is negative. A similar TS-dominated break-up can be proposed by suitable 
adjustment of the algebraic forms in (4.18) as indicated earlier but that yields a weaker and 
apparently less likely singularity. Extension of (4.22) to incorporate the vortex effect more 
fully can be made by expanding about (4.22) or solving the entire terminal system in its 
integro-differential form described above, analogous to the break-up features in (b). 

The implications of items (e), (g) in particular are discussed in the next section. 

5. REPERCUSSIONS, AND THE TYPE B INTERACTION 

The main repercussion follows from the findings in (g) in the previous section and concerns 
the next stage in the evolution of the typical vorticity tongue in the channel flow. This seems 
to occur when the spanwise scale \z\ of the streak reduces to 0(e -2 ), for the following reason. 
The relative effect in the TS wave due to its slow growth and to the influence of the induced 
vortex is of order e 4 [see Section 2] but becomes accentuated like | |[cx (T s — T) x )] typically 
as the break-up takes control. In contrast, the relative effect from classical amplitude-cubed 
contributions is smaller, of order £ 6 , but proportional to |A| 2 say, i.e. oc e 6 (T s — T) 2 , as the 
break- up is approached. Therefore the two effects become comparable when the scaled time 
T — T 3 decreases to 0(e 2 ), corresponding to the 2 -scale of 0(e -2 ) on account of (g) above. 
The next stage then, labelled Type b, has an increased vortex amplitude, with 

(u,v,w) = (u,0,0) + 0(e 2 ,e 7 ,e 5 ) (vortex) (5.1) 
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in the core, an increased TS amplitude of the form 

(«,p) = (e 2 AY,0) + 0(e 3 ,e 5 ) (TS) (5.2) 

in the two 0(e 2 )-thick viscous wall layers, and buffer zones of width 0(e) in which the skin 
friction of the vortex flow is determined. The multiple scaling also alters from that in Section 
2 to 

{d t ,d X id z ) — > (e 3 5x + , edx, e 2 dz 1 )- (5-3) 

Hence the effects of the vortex skin friction and momentum integral (oc J ) compared with 
those of the basic flow are 0(e 2 ), interacting with the new slow-temporal and warping effects 
in (5.3) and with the 0(e 2 ) amplitude-squared influence that the TS wave exerts in producing 
the vortex motion. In addition, amplitude-cubed effects of relative size 0(e 2 ) now also enter 
as a major new physical contribution in the reckoning of the TS amplitude itself. This 
stage is sketched structurally in Figure 8 and merits further study, in particular to see if the 
spanwise focusing observed in the nonlinear interaction Type a continues through the next 
stage, Type b. 

A similar inference holds if, in the spirit of item (e) in Section 4, we turn to consider 
different starting conditions. As the effective amplitude T _1 is increased the next distinct 
stage encountered is exactly that of (5.1)-(5.3). Reduced warping in the input flow, on the 
other hand, as in item (a), can lead to a much slower nonlinear-exponential evolution of the 
interaction as time increases, more similar to that found in Hall and Smith (1987). This 
slower evolution applies to the full system as well, as the footnote to item (g) mentions [see 
also additional comments on initial conditions and by-pass processes in Smith, 1988b], and 
the exponential form links up with our related work (Bennett, et ah, 1989), referred to more 
fully in the next section. 

6. NONLINEAR INTERACTION, TYPES C, D, E 

The earlier sections have shown two related Types a,b of nonlinear interaction to be 
possible. We turn now to the alternative vortex/wave interactions c,d,e which can arise from 
different input conditions as follows. 

First, we consider Types c,d, which emerge from the 3D interactive three-layer structure 
based on the e _1 length scale in x (Smith, 1980). The flow solutions in the three layers have 
the expansions 

' [t 2 U, e 5 V, e 2 W, e*P]+,y = 1 - c 2 Y+, 

[u,u,tu,p] = i [5 + e 2 Au', -eAxw,0(c 4 ),0(c 4 )],y = 0(1), 
k [e 2 U,e 5 V,t 2 W,e*P}-,y = e*Y-, 
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to the orders shown, with the unknown pressures P ± independent of Y ± in turn and now 
(x,z) — t (X, Z), t = e 3 T where again e = Re 7 . Hence the governing equations in the 
two wall layers are the 3D unsteady interactive-boundary-layer equations, 

Wx + Vy + W z ]± = 0, (6.2a) 

[U T + UU x + VU y + WU z ] ± = -P£(X, Z, T) + £/± y , (6.26) 

[Wr + UW X + VW Y + WW Z }± = -P±(X,Z,T) + Wp Y (6.2c) 

subject to no slip at Y ± = 0 and 

U* ~ y ± =F A(X, Z, T), W ± -4 0, as Y ± - oo, (6.2 d) 

with —A representing the unknown core-flow displacement. The pressure-displacement law 
necessary to close the system comes from the core-flow dynamics (c.f. Section 2), giving 

P + - P~ = I A xx (6.2c) 

where again I = 1/30. The 3D fully nonlinear formulation (6.2a-e) is the counterpart of the 
3D triple-deck formulation in boundary layers (see e.g., Hall and Smith, 1984, 1988) and in 
general it requires a computational treatment (e.g., Smith, 1989), while a possibility for its 
ultimate behavior is the finite-time singularity of Smith (1988a). 

The emergence of vortex/TS nonlinear interaction from within the system (6.2) is ad- 
dressed next, along the lines of the simultaneous work by Smith (1988b), although here 
temporal rather than spatial evolution is of concern. An order-of-magnitude argument es- 
tablishes the main possibilities. Thus, if the warping factor is 0(~ d z ) then for an input 
wave amplitude 0(6) the relative effects present, compared with 2D theory, are 

fl 2 ,h\h 2 /3 2 k-\k, (6.3 a-d) 

due to warping, to the disturbance size, to the induced vortex-shear influence, and to the 
slow temporal dependence ( k ) of the vortex, respectively. Here the third contributions stems 
from the feature that the velocity W of the induced vortex for Y of order unity is 0(6 2 /?), 
for either wall layer, and it remains at that level (apart from a logarithmic factor) in the 
vortex-dominated outer zone or buffer which is induced where Y ~ k~ 2 is large, for an 
unsteady- viscous force balance in the vortex motion. In consequence, mass conservation 
implies that the normal velocity V ~ k~\h 2 p 2 , and then the streamwise momentum balance 
indicates that U ~ k ih 2 j3 2 . This therefore yields the shear effect in (6.3c). So significant 
interaction can occur at low amplitudes (6) in either of the two cases 

/3 ~ 6, k ~ 6 2 , 6 < 1 [Type c], (6.4) 
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[Type d]. (6.5) 

Other special cases cannot be discounted but the Types c,d appear to be substantial ones. 
The former, Type c, is summarized first below. 

Type c’s double structure (with Y ~ 1 and Y ~ h~ l ) and its solution expansions proceed 
as guided by (6.4) and by comparison with Smith (1988b) and Hall and Smith (1989). The 
nonlinear interaction produced is found after some working to be governed by the induced- 
vortex equations 

V v +Wj = 0, (6.6a) 

U<+V = Uyy, ( 6 . 66 ) 

Wj + dt{\P\ 2 )Y~ 2 = Wyy (6.6c) 

subject to 

Uy -» 0, VT — > 0 as F — > oo, (6.6d) 

U ~ A 3 F, V -» 0, W — ► ^(| P\ 2 ){-£nY + <£} as F -► 0 + , (6.6e) 

coupled with the TS-wave equation 

Aj = {C! + c 3 (A+ + A3 )}A + c a ~Aj- z + c 5 A\A\ 2 (6.6/) 

for ]4(J,<). Here the coefficients in (6.6e,f) can be deduced from Sections 2-4 and Hall and 
Smith (1984), T + - ~P~ = -a 2 LA from (6.2e), and the scalings behind this interaction are 
given by 

dx — * dj + h 2 &ii dz > hdj , (6.7a) 

(P, A) = (P, A)exp(iaX — iLlT) + c.c. [wave], (6.76) 

( U,V,W ) -> (6 _1 F + hU , 6 3 V, /i 3 VT), T = [vortex], (6.7c) 

while, again as before, A 3 (z,t) denotes the correction to the skin friction. The equations 

(6.6a-e) apply at each wall (±). Lower-order logarithmic terms involved in (6.7c) have no 
influence on the nonlinear interaction (6.6a-f). The solution properties would seem to be 
analogous to those in Smith (1988b) and Section 4, implying that Type c can act as a 
fore-runner to the full interactive-boundary-layer interaction (6.2a-e) or to Type d below. 

The Type d nonlinear interaction in (6.5) is much stronger in the sense that it produces 
0(1) effects, i.e., locally the mean flow is completely altered from the original PPF form, 
despite the smallness of the input TS wave amplitude h. Here again the details are largely 
as in Smith (1988b) (see also Hall and Smith, 1989), but modified for the temporal response, 
and we find in summary the nonlinear interactive system 

Vy + W z = 0, (6.8a) 
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u- t + VUy + WUz = Uyy, (6.86) 

W t + VW Y + WW Z = Wyy, (6.8c) 

with 

Uy-+l, W -* 0, as Y -» oo, (6.8(f) 

U ~ v -* 0, -» -A; 2 ^{|P| 2 + a- 2 P z P*}, as F - 0+ (6.8c) 

for the vortex flow in the buffer zones, coupled with the TS-equation 


Pzz ~ K'KzTPz ~ <* 2 P = TGA 

for the pressures P ± and displacements oc ^A, near either wall (±). Here 


3 ^ 0 i4t' 0 

2 2Ai 0 

Ai 


/ 1«0 


G ~ {ia A v ) 3 — ~,£ 0 = — *3ft(aA v ) 3 ? 

and (6.2e) still relates and A. Also, the scalings 


(6.9a) 

(6.96) 
(6.9c, d) 


A — h 1 Y , 

(P, A) -> h'(P,A)E + c.c. [wave], (6.10a - c) 

(U, V, W) —y (h~ 1 U 1 hV,h 2 W) [vortex] 

have been applied above, with h' = h(—lnh)-*, while the wave depends on T through 
the unknown frequency fl and wavenumber a (both real) in E = exp(iaX - iSIT) but the 
induced vortex motion is independent of the fast scale T. See further comments in Smith 
(1988b). The major solution properties here are still unknown, c.f. the last-named paper, 
although the weak linear case of secondary 3D instability for a 2D input TS wave can be 
established readily, and brings in Squire-mode destabilization, as in that paper, while a 
break-up singularity at finite time is a possibility in the full nonlinear system. 

It is convenient now to reconsider the nonlinear interaction labelled Type b in Section 
5. This arises as a distinct intermediate interaction between Types a,c. Thus, if the input 
amplitude h is reduced to order e in the Type c interaction then formally, from (6.4), the 
relative time dependence k slows to 0(e 2 ), the spanwise dependence e/? reduces to 0(e 2 ) and 
the outer buffer-zone extent c 2 h 1 increases to 0(e). These suggested scales tie in exactly 
with those inferred in Section 5 from the break-up of the Type a interaction. Thus, there is 
a match between a,b and between b,c. The special character of Type b is described already 
in Section 5 and here we add only that, unlike Type c, Type b has the vortex motion being 
forced both internally and through the wall conditions by the TS amplitude. 
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The final Type, e, to be considered needs only a brief mention, as it is the short-wave/long- 
vortex nonlinear interaction of Bennett, et al. (1988), consisting of a TS wave of time scale 
0(e -3 ) interacting with a vortex of the viscous time scale 0(e -7 ). Like Type d above, Type 
e is a very strong interaction in which the mean flow itself is one of the unknown quantities 
and is substantially changed from its original form. The main point of mentioning Type e 
here is that there is a match between it and Type d, in fact, since a reduction of the input 
amplitude h to order e 2 , in Type d, increases the vortex time scale e -1 /; 1 to e -7 , in view of 
(6.5), and simultaneously increases the buffer extent y ~ e 2 h _1 to 0(1), so that the vortex 
flow occupies the entire channel. These and other features agree with the scalings for Type 
e. 

Further connections exist between the various nonlinear interactions a-e, as shown schemat- 
ically in Figure 9, e.g., a match of a with e is possible as noted in an earlier section. In 
addition, the early occurrence of each is dependent on the input amplitude ( h ) and span- 
wise dependence or warping factor (oc /?), whereas later occurrences of the interactions can 
happen either as by-passes or as the consequences of a break-up in an earlier occurring type, 
such as in Section 5 where Type b is set into action after the break-up of Type a. 

7. FURTHER COMMENTS 

As far as we are aware, the present theory is the first to show a mechanism for the 
formation of tongues, or their initiation, during channel-flow transition, starting from an 
input which is a near-2D TS wave. The mechanism involved in the Type a interaction say 
is highly interactive (see Sections 2-4), between the warped TS wave and its induced longi- 
tudinal vortex motion, and the interaction is found to be mutually reinforcing in the sense 
that the vortex flow and the warping of the TS wave amplify together very strongly. This 
is first through linear 3D secondary instability, but, beyond that, ultimately the nonlinear 
interaction leads to the vorticity-tongue formation as described in (g) in Section 4. Weaker 
warping can prevent the occurrence of such a tongue, we note, so that the present regime 
is in some respects a threshold one. The response for Type a seems quite dependent on the 
initial conditions even so and these could bring in extra, stronger, physical effects early on, 
i.e. a different structural start to the interaction could be made. One example of this is 
noted in Section 5, where Type a leads into Type b. Another is associated with the influence 
of a mean vortex motion [«( 3 \u^], or small cross-flow, present initially, which could alter 
the secondary instability properties of Section 4(b), (d) through the effects of the vortex’s 
mean skin friction and momentum integral, for instance. In general, during the current 
stage for Type a, however, Figure 7 summarizes what happens physically as the tongue is 
initiated. The vortex becomes concentrated and accentuated first at a certain spanwise lo- 
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cation, along with the enhancement of the TS amplitude and phase there, compared with 
the earlier slightly warped TS wave with weak vortex motion. In the vorticity formation, the 
integrated momentum (oc J) of the induced vortex increases rapidly in magnitude, as does 
its skin-friction factors, leading to an increasingly significant impact on the mean stream- 
wise momentum and wall shears of the total flow. Most of these theoretical features may 
tie in qualitatively at least with experimental and computational findings (see Section 1) on 
span wise focussing and eventual formation. 

Further studies indicated by the present work include the next regime in the tongue 
development (Section 5 and Figure 8), the effects of wall curvature, which may accentuate 
the vortex flow, the influence of the additional inflexional instabilities which the tongue 
behavior sets up, oblique TS waves and the application to boundary-layer transition (see 
Hall and Smith, 1988, 1989, and Smith, 1988b). The same comments apply also to Types 
b-e as well as to the 3D interactive-boundary-layer version in (6.2) which could be taken 
as a further Type, f. These interactions are of various strengths, the strongest being Types 
e,f. [Type c we note is equivalent to both cases I, II of Hall and Smith (1988c) and Type d 
corresponds to case III in the same paper.] There may be other significant types also. We 
tend to feel that it could be easy, and misleading perhaps, to miss many of the above types 
in purely computational work or crude theories. Further, it is possible that the present work 
may lead on eventually to provide insight into the structure of the so-called Lambda vortex 
and other stronger features of transitional flows. 
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Figure 1 The structure of the Type-a nonlinear interaction between a 

warped 2D TS wave and its induced longitudinal-vortex motion 
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The growth rate from (4.1 1c), when c-| r is negative. 
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Figure 4.- Spectra^ solutions computed for (4.1 a-e), (4.2): cf. Fig. 6. Time 
step AT = 0.0002, y-step Ay = 0.05, truncation at M = 40 terms 
Initial condition v ( 3 ) ( n ) = y2(i - y 2 ) > A ( n ) = ex p {-( n + 1 2)} . 
Similar features resulted for other conditions and grids. 


800 r- 



Figure 5.- Finite-difference computed solutions for (4.1a-e), (4.2). 

Here AZ = 0.1 , Ay = 0.05, AT = 0.001 . Similar trends 
were found for other conditions used. 
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Figure 6.- Spectral solutions computed for (4.1a-e), (4.2) in the special 
case c 4 = 0. The representation p = £ p( n bos(n nZ/L) 

(etc.) with L = 2 was truncated at M = 40 terms. Here 
AT = 0.0002, Ay = 0.025, Ci = 0.9017 - 0.01 189i. Initial 
condition: A( n ) = exp{-(n + I) 2 }. Graphically identical results 
were obtained for (M, AT, Ay) = (20, 0.0004, 0.05). Vortex 
velocities u( 3 ), v ( 3 ) grow similarly in amplitude. 
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Figure 7.- Sketch of the streak-like breakdown: see (g) in section 4. 
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Figure 8.- The Type-b nonlinear interaction, spanwise focussed 

from Type a as shown. Similar focussing applies to the 





Figure 9.- Schematic diagram showing where the Types a-e of nonlinear 
interaction can occur, in the normalized amplitude-spanwise 
wavenumber (h-p) plane. Type f is the triple-deck-like case (6.2) 
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